Analysis of dephasing mechanisms in a standing wave dipole trap 
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We study in detail the mechanisms causing dephasing of hyperfine coherences of cesium atoms 
confined by a far off-resonant standing wave optical dipole trap [S. Kuhr et al, Phys. Rev. Lett. 
91, 213002 (2003)]. Using Ramsey spectroscopy and spin echo techniques, we measure the reversible 
and irreversible dephasing times of the ground state coherences. We present an analytical model 
to interpret the experimental data and identify the homogeneous and inhomogeneous dephasing 
mechanisms. Our scheme to prepare and detect the atomic hyperfine state is applied at the level of 
a single atom as well as for ensembles of up to 50 atoms. 

PACS numbers: 32.80.Lg, 32.80.Pj, 42.50. Vk 



I. INTRODUCTION 

The coherent manipulation of isolated quantum sys- 
tems has received increased attention in the recent years, 
especially due to its importance in the field of quantum 
computing. A possible quantum computer relies on the 
coherent manipulation of quantum bits (qubits) , in which 
information is also encoded in the quantum phases. The 
coherence time of the quantum state superpositions is 
therefore a crucial parameter to judge the usefulness of 
a system for storage and manipulation of quantum infor- 
mation. Moreover, long coherence times are of great im- 
portance for applications in precision spectroscopy such 
as atomic clocks. 

Information cannot be lost in a closed quantum system 
since its evolution is unitary and thus reversible. How- 
ever, a quantum system can never be perfectly isolated 
from its environment. It is thus to some extent an open 
quantum system, characterized by the coupling to the en- 
vironment This coupling causes decoherence, i. e. the 
evolution of a pure quantum state into a statistical mix- 
ture of states. Decoherence constitutes the boundary be- 
tween quantum and classical physics [31 , as demonstrated 
in experiments in Paris and Boulder J^j S 13 ■ There, de- 
coherence was observed as the decay of macroscopic su- 
perposition states (Schrodinger cats) to statistical mix- 
tures. 

We can distinguish decoherence due to the progressive 
entanglement with the environment from dephasing ef- 
fects caused by classical fluctuations. This dephasing of 
quantum states of trapped particles has recently been 
studied both with ions |6l| and neutral atoms in optical 
traps 0)13 ■ In this work, we have analyzed measurements 
of the dephasing mechanisms acting on the hyperfine 
ground states of cesium atoms in a standing wave dipole 
trap. More specifically, we use the two Zeeman sublevels 
\ F = A,niF = 0) and | F = 3, niF = ) which are cou- 
pled by microwave radiation at ajhfs/27r = 9.2 GHz. 

We present our setup and the relevant experimental 



tools in Sec.ni with special regard to the coherent manip- 
ulation of single neutral atoms. Our formalism and the 
notation of the dephasing/decoherence times are briefly 
introduced in Sec. IIIII Finally, in Sees. IIVI and El we 
experimentally and theoretically analyze the inhomoge- 
neous and homogeneous dephasing effects. 



II. EXPERIMENTAL TOOLS 
A. Setup 

We trap and manipulate cesium atoms in a red de- 
tuned standing wave dipole trap. Our trap is formed 
of two counterpropagating Gaussian laser beams with 
waist 2wo = 40 pm and a power of max. 2 W per 
beam (see Fig. derived from a single Nd:YAG laser 
(A ~ 1064 nm). Typical trap depths are on the order of 
Uq = 1 mK. The laser beams have parallel linear polar- 
ization and thus produce a standing wave interference 
pattern. Both laser beams are sent through acousto- 
optic modulators (AOMs), to mutually detune them for 
the realization of a moving standing wave. This "optical 
co nvey or-belt" was introduced in previous experiments 
0, llOj and has been used for the demonstration of quan- 
tum state transportation prj . For the experiments in 
this paper, however, we do not transport the trapped 



MOT laser 
beams 



MOT 



Vacuum 
chamber 



Dipole trap 
laser beam 



Dipole trap 
laser beam 




'Electronic address: kuhr @lkb.ens.frl 



APD 

FIG. 1: Experimental setup. 
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atoms. To eliminate any heating effect arising due to 
the phase noise of the AOM drivers |1QL I12| , we used the 
non-deflected beams (O**^ order of the AOMs) to form 
the dipole trap. The AOMs are only used to vary the 
dipole trap laser intensity by removing power from the 
trap laser beams. 

Cold atoms are loaded into the dipole trap from a 
high gradient magneto-optical trap (MOT). The high 
field gradient of the MOT {dB/dz = 340 G/cm) is pro- 
duced by water cooled magnetic coils, placed at a dis- 
tance of 2 cm away from the trap. The magnetic field 
can be switched to zero within 60 ms (limited by eddy 
currents in the conducting materials surrounding the vac- 
uum chamber) and it can be switched back on within 
30 ms. Our vacuum chamber consists of a glass cell, 
with the cesium reservoir being separated from the main 
chamber by a valve. Cesium atoms are loaded into the 
MOT at random from the background gas vapor. To 
speed up the loading process, we temporarily lower the 
magnetic field gradient to dB/dz — 25 G/cm during a 
time iiow The low field gradient results in a larger cap- 
ture cross section which significantly increases the load- 
ing rate. Then, the field gradient is returned to its initial 
value, confining the trapped atoms at the center of the 
MOT. Varying tiow enables us to select a specific average 
atom number ranging from 1 to 50. The required time 
depends on the cesium partial pressure, which was kept 
at a level such that we load typically 50 atoms within 
^low = 100 ms in these experiments. 

In order to transfer cold atoms from the MOT into 
the dipole trap, both traps are simultaneously operated 
for some tens of milliseconds before we switch off the 
MOT. After an experiment in the dipole trap the atoms 
are transferred back into the MOT by the reverse proce- 
dure. All our measurements rely on counting the number 
of atoms in the MOT before and after any intermediate 
experiment in the dipole trap. For this purpose we col- 
lect their fiuorescence light by a homc-built diffraction- 
limited objective [T^ and detect the photons with an 
avalanche photodiode (APD). 

Three diode lasers are employed in this experiment 
which are set up in Littrow configuration and locked by 
polarization spectroscopy. The MOT cooling laser is sta- 
bilized to the F = 4— >F' = 3/i^' = 5 crossover transi- 
tion and shifted by an AOM to the red side of the cooling 
transition F ~ A ^ F' = b. The MOT repumping laser is 
locked to the F — 2> ^ F' — A transition, it is 7r-polarized 
and is shined in along the dipole trap axis. To optically 
pump the atoms into the \ F = A^mp = 0) state, we use 
the unshifted MOT cooling laser, which is only detuned 
by -1-25 MHz from the required F ^ A ^ F' = A tran- 
sition. This small detuning is partly compensated for 
by the light shift of the dipole trap. We shine in the 
laser along the dipole trap axis with 7r-polarization to- 
gether with the MOT repumper. We found that 80% of 
the atoms are pumped into the | F — A,mp — 0) state, 
presumably limited by polarization imperfections of the 
optical pumping lasers. 



For the state selective detection (see below) we use 
a "push-out" laser, resonant to the = 4 — > i^' = 5 
transition. It is (T+-polarized and shined in perpendicular 
to the trapping beams (z-axis in Fig. 

To generate microwave pulses at the frequency of 
9.2 GHz we use a synthesizer (Agilent 83751 A), which is 
locked to an external rubidium frequency standard (Stan- 
ford Research Systems, PRSIO). The amplified signal 
{P — +36 dBm — 4.0 W) is radiated by a dipole an- 
tenna, placed at a distance 5 cm away from the MOT. 

Compensation of the earth's magnetic field and stray 
fields created by magnetized objects close to the vac- 
uum cell is achieved with three orthogonal pairs of coils. 
For the compensation, we minimize the Zeeman split- 
ting of the hyperfine ground state mp manifold which is 
probed by microwave spectroscopy. Using this method 
we achieve residual fields of B^cs < 0.4 jiT (4 mC). The 
coils of the z-axis also serve to produce a guiding field, 
which defines the quantization axis. 



B. State selective detection of a single neutral 
atom 

Sensitive experimental methods had to be developed in 
order to prepare and to detect the atomic hyperfine state 
at the level of a single atom. State selective detection is 
performed by a laser which is resonant with the F — 
4 — > F' = 5 transition and thus pushes the atom out of 
the dipole trap if and only if it is in = 4. An atom 
in the _F = 3 state, however, is not influenced by this 
laser. Thus, it can be transferred back into the MOT 
and be detected there. Although this method appears 
complicated at first, it is universal, since it works with 
many atoms as well as with a single one. Other methods, 
such as detecting fiuorescence photons in the dipole trap 
by illuminating the atom with a laser resonant to the 
F = A ^ F' = 5 transition, failed in our case because 
the number of photons detected before the atom decays 
into the F = 3 state is not sufficient. 

In order to achieve a high efficiency of the state selec- 
tive detection process, it is essential to remove the atom 
out of the dipole trap before it is off-resonantly excited 
to F' = 4 and spontaneously decays into the _F = 3 state. 
For this purpose, we use a cr"'"-polarized push-out laser, 
such that the atom is optically pumped into the cycling 
transition \ F = A,mp = A) ^ \ F' = 5, mp — 5). In our 
setup, the polarization is not perfectly circular, since for 
technical reasons we had to shine in the laser beam at an 
angle of 2° with respect to the magnetic field axis. This 
entails an increased probability of exciting the F' — A 
level from where the atom can decay into the F = 3 
ground state. To prevent this, we remove the atom from 
the trap sufficiently fast by shining in the push-out laser 
from the radial direction with high intensity (///o ~ 100, 
with wq — 100 pm, P = 30 pW, where /g — 1-1 mW/cm^ 
is the saturation intensity). In this regime its radiation 
pressure force is stronger than the dipole force in the ra- 
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FIG. 2: State selective detection of a single atom. The graphs 
show the fluorescence signal of the atom during state prepa- 
ration and detection, binned in time intervals of 5 ms. The 
bars above the graphs show the timing of the lasers. Graphs 
(a)(i) and (b)(i) show the signals of a single atom, prepared 
in F = 3 and F — 4, respectively. Graphs (a)(ii) and (b)(ii) 
show the added signal of about 150 events. 



dial direction, such that we push out the atom within half 
the radial oscillation period 1 ms). In this case, the 
atom receives a momentum corresponding to the sum of 
all individual photon momenta. This procedure is more 
efficient than heating an atom out of the trap, which oc- 
curs when the radiation pressure force of the push-out 
laser is weaker than the dipole force, and the atom per- 
forms a random walk in momentum space while absorb- 
ing and emitting photons. 

If we adiabatically lower the trap to typically 0.12 mK 
prior to the application of the push-out laser, we need on 
average only 35 photons to push the atom out of the trap. 
This number is small enough to prevent off-resonant ex- 
citation to F' = 4 and spontaneous decay to F = 3. 
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FIG. 3: Atom counting. Initial and final numbers of atoms 
are inferred from their fluorescence in the MOT. Shown are 
the integrated APD counts binned in time intervals of 1 ms 
and accumulated over 10 repetitions with 20 atoms each. 



A typical experimental sequence to test the state se- 
lective detection is shown in Fig. El First, the atom is 
transferred from the MOT into the optical dipole trap. 
Using the cooling and repumping laser of the MOT, we 
optically pump the atom either in the F = 3 (Fig.|5^) or 
the F ~ 4 hyperfine state (Fig. I^b). The push-out laser 
then removes all atoms in F = 4 from the trap. Any 
remaining atom in F = 3 is transferred back into the 
MOT, where it is detected. 

Figures 13^(1) and |2t>(i) show the signals of a single 
atom, prepared in F = 3 and F = 4, respectively. Our 
signal-to-noise ratio enables us to unambiguously detect 
the surviving atom in F = 3, demonstrating the state- 
selective detection at the single atom level. We performed 
157 repetitions with a single atom prepared in F = 3 
and found that in 153 of the cases the atom remains 
trapped, yielding a detection probability of 97.5^2 o*^. 
Similarly, only 2 out of 167 atoms prepared in F = 4 
remain trapped, yielding 1.2^q g%. The asymmetric er- 
rors are the Clopper-Pearson 68% confidence limits [l^ . 
These survival probabilities can also be inferred by di- 
rectly adding the signals of the individual repetitions and 
by comparing the initial and final fluorescence levels in 
the MOT, see Figs.|li(ii) and|l)(ii). 

All following experiments are performed in the 
same way. We initially prepare the atoms in the 
|F = 4, mi? = 0) state and measure the population 
transfer to \ F — S^mp = 0) induced by the microwave 
radiation. After the application of one or a sequence of 
microwave pulses, the atom is in general in a superposi- 
tion of both hyperfine states, 

\^) = C3 \ F = 3,mF ^ 0) + C4 \ F = 4,mF ^ 0) , (1) 

with complex probability amplitudes C3 and C4. Our de- 
tection scheme only allows us to measure the population 
of the hyperfine state F = 3: 



„ 1,9 W + 1 



(2) 



where w is the third component of the Bloch vector, see 
below. The number P3 is determined from the number of 
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FIG. 4: Rabi oscillations on the \ F — i, mp ~ 0) 
I F = 3, mp = ) clock transition recorded at a trap depth 
Uo ~ 1.0 mK. Each data point results from 100 shots with 
about 60 initial atoms. The line is a fit according to Eq. Q. 



atoms before (iVinitiai) and after (iVfi„ai) any experimen- 
tal procedure in the dipole trap. A^initiai and iVfinai are 
inferred from the measured photon count rates, Cinitiai, 
Cfinai and Cbackgr (see Fig. El : 
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The fluorescence rate of a single atom, Ciatom, is mea- 
sured independently. From the atom numbers we obtain 
the fraction of atoms transferred to F = 3, 



final 



(5) 



The measured number of atoms, iVjnitai can be larger than 
the actual number of atoms in the dipole trap, since we 
lose atoms during the transfer from the MOT into the 
dipole trap (see below). 



C. Rabi oscillations 

We induce Rabi oscillations by a single resonant mi- 
crowave pulse at the maximum RF power. For the graph 
of Fig. ^ we varied the pulse length from ps to 225 ps 
in steps of 5 ps. Each point in the graph results from 100 
shots with about 60 ± 10 atoms each. The corresponding 
statistical error is below 1% and is thus not shown in the 
graph. The error of the data points in Fig. 0] is domi- 
nated by systematic drifts of the storage probability and 
efficiencies of the state preparation and detection. Since 
'w{t) = — cos ^Rt, we fit the graph with 



C 



1 — cos riRt) , 



(6) 



which yields a Rabi frequency f7r{/27r ~ (14.60 ± 
0.02) kHz. Note that this Rabi frequency is higher than 
the one used later in this report (10 kHz) because we 



changed the position of the microwave antenna for practi- 
cal reasons. The maximum population detected in F = 3 
is C = (60.4±0.7)%. This reduction from 100% is caused 
by two effects. First, when we use many (> 40) atoms 
at a time, up to 20% of the atoms are lost during the 
transfer from the MOT into the dipole trap due to inelas- 
tic collisions, as verified in an independent measurement. 
The remaining losses arise due to the non-perfect optical 
pumping process. 



HI. PHENOMENOLOGICAL DESCRIPTION OF 
DECOHERENCE AND DEPHASING 

In our experiment we observe quantum states in an 
ensemble average, and decoherence manifests as a decay 
or dephasing of the induced magnetic dipole moments. 
It is useful to distinguish between homogeneous and in- 
homogeneous effects. Whereas homogeneous dephasing 
mechanisms affect each atom in the same way, inhomoge- 
neous dephasing only appears when observing an ensem- 
ble of many atoms possessing slightly different resonance 
frequencies. As we will see later, the most important 
difference between the two mechanisms is the fact that 
inhomogeneous dephasing can be reversed, in contrast to 
the irreversible homogeneous dephasing. 

The interaction between the oscillating magnetic field 
component of the microwave radiation, B cosujt, and the 
magnetic dipole moment, fi, of the atom is well approxi- 
mated by the optical Bloch equations J/Tj: 



n X u 



(7) 



with the torque vector ft = (JIr, 0, S) and the Bloch 
vector u = {u,v,w). Here, JIr = fsB/Ti is the Rabi 
frequency and 5 ^ lo — loq is the detuning of the mi- 
crowave from the atomic transition frequency uj^. In 
the following, the initial quantum state \ F — A^mp = 
corresponds to the Bloch vector u = (0,0,-1), whereas 
I F = 3,mp = 0) corresponds to it = (0,0, 1). 

We include the decay rates as damping terms into the 
Bloch equations and use a notation of the different times 
for population and polarization decay similar to the one 
of nuclear magnetic resonance: 



= 5{v)~ 

= -fiR (V 



T2 
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(w) - Wst 



Ti 



(8a) 
(8b) 
(8c) 



where (...) denotes the ensemble average. The total ho- 
mogeneous transverse decay time T2 is given by the po- 
larization decay time T2' and the reversible dephasing 



time T2 



1 



1 

TV 



(9) 
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Inhomogeneous dephasing (Tj ) occurs because the atoms 
may have different resonance frequencies depending on 
their environment. Thus the Bloch vectors of the in- 
dividual atoms precess with different angular velocities 
and lose their phase relationship, they dephase. In our 
case, inhomogeneous dephasing arises due to the energy 
distribution of the atoms in the trap. This results in a 
corresponding distribution of light shifts because hot and 
cold atoms experience different average trapping laser in- 
tensities. 

The longitudinal relaxation time, Ti, describes the 
population decay to a stationary value Wst- In our case, 
Ti is governed by the scattering of photons from the 
dipole trap laser, which couples the two hyperfine ground 
states via a two-photon Raman transition. This effect is 
suppressed due to a destructive interference effect yield- 
ing relaxation times of several seconds (see Sec. lVB| . We 
do not include losses of atoms from the trap in the decay 
constants, which occur on the same timescale. 



IV. INHOMOGENEOUS DEPHASING 

We measure the transverse decay time T2 by perform- 
ing Ramsey spectroscopy, which consists of the applica- 
tion of two coherent rectangular microwave pulses, sepa- 
rated by a time interval t The initial Bloch vector 
Mo = (0,0,-1) corresponds to an atom prepared in the 
\F = A, rap — 0)-state. A 7r/2-pulse rotates the Bloch 
vector into the state (0,-1,0), where the atom is in a 
superposition of both hyperfine states. The Bloch vector 
freely precesses in the itw-plane with an angular frequency 
5. Note that i5 has to be small compared to the Rabi fre- 
quency and the spectral pulse width, such that the pulse 
can be approximated as near resonant, and complete pop- 
ulation transfer can occur. After a free precession during 
t, a second 7r/2-pulse is applied. The measurement of the 
quantum state finally projects the Bloch vector onto the 
w axis. 

We recorded Ramsey fringes for two different dipole 
trap depths, C/o = 0.1 niK and 0.04 mK (see Fig.O. Each 
point in the graph of Fig. corresponds to 30 shots with 
about 50 trapped atoms per shot, yielding errors (not 
shown) of less than 1%. The quoted values for [/q are 
calculated from the measured power and the waist of the 
dipole trap laser beam and have an estimated uncertainty 
of up to 50%. We initially transfer the atoms from the 
MOT into a deeper trap (C/q > 1 mK) to achieve a high 
transfer efficiency. When the MOT is switched off, we 
adiabatically lower the trap depth using the AOMs. 

Our Ramsey fringes show a characteristic decay, which 
is not exponential. This decay is due to inhomogeneous 
dephasing, which occurs because after the first 7r/2- 
pulse, the atomic pseudo-spins precess with different 
angular frequencies. In the following, we derive analytic 
expressions for the observed Ramsey signal and we show 
that the envelope of the graphs of Fig. [3 is simply the 
Fourier transform of the atomic energy distribution. 
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FIG. 5: Ramsey fringes recorded for two different trap depths 
(a) Uo = 0.1 mK and (b) 0.04 mK. Their decay with time 
constants T2 = 4.4±0.1 ms and 20.4±1.1 ms, respectively, is 
governed by inhomogeneous dephasing caused by the energy 
distribution in the trap. Each data point results from 30 
shots with about 50 initial atoms. The damped oscillation 
is a fit with P3, Ramsey (t) and the envelopes are the functions 
B±Aa{t,T2) (see Eqs. ll^l. 



A. Differential light shift and decay of Ramsey 
fringes 

The light shift of the ground state due to the Nd:YAG 
laser is simply the trapping potential 



(10) 



The detuning of the Nd:YAG-laser from the D-line of 
an atom in F = 4 is 9.2 GHz less than for an atom in 
F = 3. As a consequence, the F — A level experiences 
a slightly stronger light shift, resulting in a shift of the 
F — Z ^ F = A microwave transition towards smaller 
resonance frequencies. This differential light shift, (5o, 
can be approximated as 



Mo = ?7o(Acff) - C/o(Aoff + t^iits), 



(11) 



where Aoff = -1.2x10"^ F is an effective detuning, taking 
into account the weighted contributions of the Di and D2 
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lines Whfs = 2.0 X 10^ r is the ground state hyperfine 
splitting. Since Whfs ^eS, we find that the differential 
light shift is proportional to the total light shift Uq, 



(12) 



with a scaling factor ri = Whfs/Acff = 1.45 X 10"'*. For 
atoms trapped in the bottom of a potential of J7o = 1 mK, 
the differential light shift is 6q = —2tt x 3.0 kHz. 

In the semiclassical limit, i. e. neglecting the quantized 
motion of the atom in the dipole trap potential, the free 
precession phase accumulated by an atomic superposition 
state between the two 7r/2-pulses depends on the average 
differential light shift only. In the following, we calcu- 
late the expected Ramsey signal using this semiclassical 
approach and obtain simple analytical expressions. Fur- 
thermore, we verified the validity of the presented model 
by performing a quantum mechanical density matrix cal- 
culation (not presented here) which agrees to within one 
percent with the semiclassical results. The small devia- 
tion can be attributed to the occurrence of small oscil- 
lator quantum numbers riosc — 5 in the stiff direction of 
the trap. 

Note that, strictly speaking, our model of a time- 
averaged differential light shift is only correct if the atom 
carries out an integer number of oscillation periods in 
the trap between the two 7r/2- pulses. However, we have 
checked that the variable phase accumulated during the 
remaining fraction of an oscillation period does not cause 
a measurable reduction of the Ramsey fringe contrast and 
can therefore be neglected. 

Since a hot atom experiences a lower laser intensity 
than a cold one, its averaged differential light shift is 
smaller. The energy distribution of the atoms in the 
dipole trap obeys a three-dimensional Boltzmann distri- 
bution with probability density |l^ 0| 



piE) = 



exp 



E 



(13) 



Here E = iJkin -I- J7 is the sum of kinetic and potential 
energy. In a harmonic trap the virial theorem states that 
the average potential energy is half the total energy, U — 
E/2. Thus, the average differential light shift for an atom 
with energy E is given by: 



^.(i^)=^o + f 



(14) 



where (5o < is the maximum differential light shift. As a 
consequence, the energy distribution p{E) yields, except 
for a factor and an offset, an identical distribution a((5is) 
of differential light shifts Q : 

a{5is) = — (<5is - <^o)' exp [~K{5^, - So)] (15) 



with 



K 



2h 



(16) 



Note that this distribution is only valid in the regime 
k^T Uo, since the virial theorem was applied for the 
case of a harmonic potential. 

To model the action of the Ramsey pulse sequence, 
we express the solutions of Eq. ^ as rotation matrices 
acting on the Bloch vector. The Ramsey sequence then 
reads 

MRamsoy(i) = ' *frco(5, ' 07r/2 ' Uq, (17) 

with the matrices describing the action of a 7r/2-pulse, 

'l O' 
1 
,0 -1 



6^/2 =10 1 



(18) 



and the free precession around the w-axis with angular 
frequency 5 during a time interval i. 



*frco('5,i) 



cos0((5, i) sin 0(5, t) 
— sin0((5, i) cos(/)((5, t) 



(19) 











The total precession angle, 4>{6,t), represents the accu- 
mulated phase during the free evolution of the Bloch vec- 
tor, 0(t) = 6{t')dt' . The detuning S{t) may in general 
vary spatially and in time, depending on the energy shifts 
of the atomic levels. 

If the Bloch vector is initially in the state Uq — 
(0, 0,-1), we obtain from Eq. (fTTjl : 



WRamsoy(i) = COS St, 



(20) 



where S ~ lu — loq is the detuning of the microwave ra- 
diation with frequency uj from the atomic resonance luq. 
In order to see the Ramsey fringes, we purposely shift 
LO with respect to the ground state hyperfine splitting, 
^hfsi by a small detuning, Jgynth, set at the frequency 
synthesizer 



(21) 



The atomic resonance frequency ujo is modified due to 
external perturbations 



^0 = ^^hfs + Sis + Sb, 



(22) 



where S\s is the energy dependent differential light shift, 
Sb is the quadratic Zeeman shift. 

Now, the inhomogeneously broadened Ramsey signal 
is obtained by averaging over all differential light shifts 

S\s- 

/•oo 

WRamsoy,mh(0 = J a{Sls) 

X COS [((5synth - Sb ~ Sls)t] dSls- 

Eq. (|23|) shows that the shape of the Ramsey fringes is the 
Fourier(-Cosine)-Transform of the atomic energy distri- 
bution. Note that in the above integral we have set the 
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upper integration limit to oo, instead of the maximum 
physically reasonable value, 5q/2, to obtain the analytic 
solution 

^Ramsey ,inh 

with the sum of the detunings, 

6' = (^synth - - <^0 (25) 

and a time- dependent amplitude a(i, T2 ) and phase shift 



a{t,T;) = 



1 + 0.95 



^2 



-3/2 



and 



K(i,r*) = -Sarctan ( 0.97— 

J 9 



(26) 



(27) 



Despite this non-exponential decay, we have intro- 
duced the inhomogeneous or reversible dephasing time 



T2 as the 1/e-time of the amplitude a(t): 



T; = Ve2/3 - 1 K = 0.97 



2n 

rfkBT 



(28) 



Thus, the reversible dephasing time T2 is inversely pro- 
portional to the temperature of the atoms. 

The phase shift K(i,T2*) arises due to the asymmetry 
of the probability distribution q;((5is). The hot atoms in 
the tail of the energy distribution dephase faster than 
the cold atoms, due to their larger spread. The fact that 
these hot atoms no longer contribute to the Ramsey sig- 
nal results in a weighting of the mean S\s towards larger 
negative values. 

To fit our experimental data, we derive the following 
expression from Eq. 124|l . 



-Oi, Ramsey 

X A cos [S't + K{t, T2) + (p] , 



(29) 



where the amplitude A and the offset B account for the 
imperfections of state preparation and detection. The 
other fit parameters are S', Tj*, and a phase offset tp. 

The corresponding fits are shown in Fig. [S] and the re- 
sulting fit parameters are summarized in TableQ] For the 
two graphs, the maximum population detected in F — 
-Pa, max = A + B is only about 60% and 30%, respectively. 
The reduction to 60% in Fig. |Sfa) is again due to imper- 
fections in the optical pumping process and due to losses 
by inelastic collisions, as discussed in Sec. Ill CI The ad- 
ditional reduction in Fig.|3Jb) occurs during the lowering 
of the trap to Uq = 0.04 mK, where another 50% of the 
atoms are lost. Note, however, that the fringe visibility 



V 



A 
B 





Fig. Ha) 


Fig.Hb) 


Uo (est.) 


0.1 mK 


0.04 mK 


<5synth/27r 


2250 Hz 


1050 Hz 


A 


28.7 ± 0.5 % 


13.6 ± 0.1 % 


B 


30.5 ± 0.1 % 


13.8 ± 0.1 % 


S'/2tv 


2133.7 ± 1.5 Hz 


722.5 ± 0.5 Hz 




0.35 ±0.02 


0.13 ±0.03 




4.4 ± 0.1 ms 


20.4 ± 0.6 ms 



TABLE I: Fit parameters extracted from the Ramsey fringes 
of Fig. El using Eq. 



is not impaired by these imperfections. From the fit pa- 
rameters we obtain V = 0.97 ± 0.01 and V = I.OOIq 03 
for the two cases. 

As a check of consistency, we calculate the differential 
light shift (5o from the fitted detuning 5' and the experi- 
mental values of 5^ and (Jsynth, 



^0 — <5syntli — & — (5'. 



(31) 



(30) 



The calculated quadratic Zeeman shift in the externally 
applied guiding field of B = 97.9 ± 1.5 yT is 5b/2tt = 
412 ± 13 Hz, where the error is due to the uncertainty 
of the calibration. We obtain 5q/2tt — —268 ± 13 Hz 
and (5o/27r = —78 ± 13 Hz. From the values of we 
can formally deduce the potential depth corresponding 
to Uo = 0.090 ± 0.004 mK and Uq = 0.026 ± 0.004 mK, 
which almost match the expected trap depths estimated 
from the dipole trap laser power assuming purely linear 
polarization. The discrepancy for the lowest trap depth 
could arise from the fact that the energy distribution is 
truncated at E = Uq, since we have lost the atoms with 
the highest energy during the lowering of the trap. This 
truncation will reduce the effective and thus yield a 
smaller trap depth. 

Finally, the phase offset tp occurs because the Bloch 
vector precesses around the w-axis even during the ap- 
plication of the two 7r/2-pulses. In contrast, our ansatz 
of Eq. H17fl takes into account only the free precession in 
between the two pulses. The additional precession angle 
amounts to 1^9 = 21^^/26'. Given = 16 ys and the 
fitted value of S' we obtain tp = 0.42 for Fig. [SJa) and 
ip = 0.14 for Fig.|3fb), which is close to the fitted values 
of Table n 



V. HOMOGENEOUS DEPHASING 
MECHANISMS 

A. Spin echoes 

The inhomogeneous dephasing can be reversed using 
a spin echo technique, i. e. by applying an additional 
TT-pulse between the two Ramsey 7r/2-pulses. Although 
originally invented in the field of nuclear magnetic reso- 
nance |l3, this technique was recently also employed in 
optical dipole traps 
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FIG. 6: Spin echoes. Shown are spin echoes recorded for three different trap depths, (a) Uo = 1.0 mK, (b) 0.1 mK and (c) 
0.04 mK. We observe a decrease of the maximum spin echo amphtude with increasing time of the vr-pulse, with longer decay 
times in lower trap depths. All spin echoes are fitted using Eq. (I37I I. In (a) and (b), the first curve is a Ramsey signal, recorded 
with otherwise identical parameters. 



We recorded echo signals in three different trap depths, 
Uq = 1.0 mK, 0.1 mK and 0.04 mK for different times of 
the TT-pulse, r^r (see Fig.El). We observe that the visibility 
of the echo signals decreases if we increase r^r. A slower 
decrease of the visibility is obtained in lower traps. For 
Uo — 0.04 mK, — 200 ms, we even observed oscilla- 
tions that reappear at t = 400 ms. 

In order to interpret these results, we first model the 
action of the microwave pulses for the spin echo, similar 
to the discussion in Sec. IIVI After the first 7r/2-pulse 
at t = 0, all Bloch vectors start at tt(0) — (0,-1,0). 
Due to inhomogeneous dephasing, the Bloch vectors ro- 
tate at slightly different frequencies around the w-axis. 
A TT-pulse at time r^r rotates the ensemble of Bloch vec- 
tors around the u axis by 180° and induces a complete 
rephasing at 2T7r in the state u{2tt^) — (0,1,0). The 
corresponding matrix equation reads 

Moclio(i) = 07r/2 • 'J'frco('5, t - T^) ■ Q^, ■ 

•*froo(<5, r^) • 6^/2 ■ Mo, (32) 

where we defined 

/l 0\ 

^ -1 . (33) 



Here, r^r is the time between the first tt/2- and the tt- 
pulse, and t > t.^ is the time of the second tt/2 pulse. As 
a result of Eq. (|32|l , we obtain 

Wccho{t) = - cos[S{t - 2r„)]. (34) 

We calculate the shape of the inhomogeneously broad- 
ened echo signal, u'ccho,inh(^); by integrating over all dif- 
ferential light shifts 6\s 

/•oo 

Wccho,inh(0 = - / &{^\s) X 
J So 

X COS [((5synth - Sis - <5B)(i - 2r^)] dSis- (35) 

The integration yields a result similar to Eq. (|23|l . 

^«ccho,inii(i) = -a{t - 2t^) 

xcos[(5'(t-2r^) + K(t-2T^)], (36) 

with amplitude a{t) and phase shift K{t) as defined in 
Eqs. and ^7^. Eq. shows that the amplitude of 
the echo signal regains its maximum at time 2T7r. Finally, 
the population in F = 3 reads: 

P3,echo(i) = B-a{t~2T^,T;) 

xAcos [5'{t- 2t^) + k(< - 2t^, T^) -f t/-]. (37) 
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Fig. Ha) 


Fig.Hb) 


Fig. He) 


Uo (est.) 

T2' 
Ti (calc.) 


1.0 mK 
0.86 ±0.05 ms 
10.2 ± 0.4 ms 
8.6 s 


0.1 iTiK 
2.9 ±0.1 ms 
33.9 ±1.0 ms 

86 s 


0.04 mK 
18.9 ± 1.7 ms 
146.2 ±6.6 ms 
220 s 



TABLE II: Summary of dephasing times. T2 and T2' are 
obtained from the echo signals of Fig. |S| Ti = Fp^Jj^j^j^ is 
calculated using Eq. 16311 . 



This equation is used to extract dephasing times from 
all spin echoes of Fig. El The average values are listed in 
Table ^ where Tj* was obtained by averaging over the 
respective datasets. From the amplitude, A, and offset, 
B, of each echo signal we calculate the visibility V — 
A/B, plotted in Fig. [7| as a function of r^r. The phase 
shift ip accounts for slow systematic phase drifts during 
the spin echo sequence. 

So far, we considered the detuning as constant during 
the experimental sequence. We now include in our model 
a time- varying detuning, S{t), in order to account for a 
stochastic variation of the precession angles of the Bloch 
vector, 



5{t) dt and 



2t, 



5{t) dt, (38) 



before and after the 7r-pulse. The phase difference 4)2 — 4>i 
is expressed as a mean difference of the detuning. 



A(5 



(39) 



The Bloch vector at time 2tt^, when the inhomogeneous 
dephasing has been fully reversed, reads 

itccho(A(5, 2t^) = • ^>frcc(<5 + A5, r^) • e^r ■ 

•«'ft-co(^,r^) • 6^2 • Wo, (40) 



which results in 

Wocho(A(5, 2t^) 



- cos(A(5 T^). 



(41) 



For a Gaussian distribution of fluctuations with mean 
A(5 = and variance (t(t^)^, 



Pr,(A5) 



1 



cr(T^)V27r 



exp 



(A5)2 



(42) 



the average w-component of the Bloch vector is calcu- 
lated, 



a) 
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FIG. 7: Decay of the spin echo visibility, extracted from the 
signals of Fig. |S] The fits (red lines) are the Gaussians of 
Eq. 1451 . The dashed and dotted lines are the best and worst 
case predictions inferred from the measured pointing instabil- 
ity of the trapping laser shown in Fig. |5] 



For comparison with the experimental values, we fit 
the spin echo visibility of Fig. [7| with a Gaussian, 



y(2r,) =Coexp 



^2^2 

— T (7 

2 TT CXp 



(45) 



with a time- independent detuning fluctuation ficxp- We 
define the homogeneous dephasing time T2 as the 1/e 
decay time of the spin echo visibility: 



v{2T2') = a 



be 



^ T2' = 



-'exp 



(46) 



'^ccho , horn ( 2 Tti- ) 



-cos(A(5t^) Pt,(A(5) dA(5 



exp 



1-^2 



Thus, the spin-echo visibility, V , yields 

F(2t,) = T/o exp [-irXT,)2] . 



(43) 
(44) 



B. Origins of irreversible dephasing 

Candidates for irreversible dephasing mechanisms 
include intensity fluctuations (1) and pointing instability 
of the dipole trap laser (2), heating of the atoms 
(3), fluctuating magnetic fields (4), fluctuations of 
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the microwave power and pulse duration (5) and spin 
relaxation due to spontaneous Raman scattering from 
the dipole trap laser (6). 





Uo 


1.0 mK 


0.1 mK 


0.04 mK 


(Toxp (meas.) 


22.0 ± 0.9 Hz 


6.6 ± 0.2 Hz 


1.54 ± 0.07 Hz 


(1) 


intensity 
fluctuations 


5.9 Hz 


0.67 Hz 


0.17 Hz 


(2) 


pointing instability 
best case 
worst case 


10.6 Hz 
21.6 Hz 


2.4 Hz 
6.7 Hz 


1.3 Hz 
3.7 Hz 


(3a) 
(3b) 


heating 

(upper limit) 

photon scattering 
rTp(T^)/27r 


5.3 Hz 
4.5 Hz 


1.6 Hz 
1.5 Hz 


2.0 Hz 
1.4 Hz 


(4) 


magnetic field 
fluctuations 
<Tb(T2')/27r 


1.7 Hz 


0.35 Hz 


0.17 Hz 













TABLE III: Summary of dephasing mechanisms. Shown are 
the fluctuation amplitudes <j{T2)/2n. 



(1) Intensity fluctuations of the trapping laser. The 
intensity fluctuations are measured by shining the laser 
onto a photodiode and recording the resulting voltage as 
a function of time. From this signal we calculate cr(r^)^ 
by means of the Allan variance, defined as j23| : 



{Xr,k 



k+1 — Xr,k) 



(47) 



k=l 



Here Xr^k denotes the average of the photodiode voltages 
over the fc-th time interval r, normalized to the mean 
voltage of the entire dataset. The resulting Allan devi- 
ation (7a is a dimensionless number which expresses the 
relative fluctuations. They directly translate into fluctu- 
ations cr(r) of the detuning, 



cr(T) = \/28a(7A[T). 



(48) 



The factor of %/2 arises because (T(r) is the standard 
deviation of the difference of two detunings with stan- 
dard deviation crA('7) each. The maximum differential 
light shift (5o in Eq. (|48|l is calculated according to 
Eq. (|31|l using the measured values of 5' . As a result 
we find relative intensity fluctuations of crA(T) < 0.2% 
(see Fig. ISJ. The corresponding absolute fluctuation 
amplitudes cr(T2')/27r (shown in Table HIHl are too weak 
to account for the observed decay of the spin echo 
visibility. 

(2) Pointing instability of the trapping laser. Any 
change of the relative position of the two interfering laser 
beams also changes the interference contrast, and hence 
the light shift Sq. These position shifts can arise due to 
shifts of the laser beam itself, due to variations of the 
optical paths e. g. from acoustic vibrations of the mirrors 
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o 
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FIG. 8: Allan deviation of the intensity fluctuations accord- 
ing to Eq. JTfl . 
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FIG. 9: Measuring the pointing instability, (a) The dipole 
trap beams having a frequency difference of l\v = 10 MHz 
are overlapped on a fast photodiode. (b) Allan deviation of 
the amplitude of the resulting beat signal. 



or from air flow. In order to measure the pointing in- 
stabilities we mutually detune the two dipole trap beams 
by Av = 10 MHz using the AOMs and overlap them on 
a fast photodiode (see Fig. (SJa)). The amplitude of the 
resulting beat signal directly measures the interference 
contrast of the two beams and is thus proportional to the 
depth of the potential wells of the standing wave dipole 
trap. We used a network analyzer (HP 3589A) operated 
in "zero span" mode to record the temporal variation of 
the beat signal amplitude within a filter bandwidth of 
10 kHz. 

The resulting Allan deviation of the beat signal ampli- 
tudes is shown in Fig. |^b). The lower curve shows the 
signal in the case of well overlapped beams, whereas for 
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the upper curve, we purposely misaligned the beams so 
that the beat signal amplitude is reduced by a factor of 
2. In the latter case variations of the relative beam posi- 
tion cause a larger variation of the beat signal amplitude, 
since the beams overlap on the slopes of the Gaussian 
profile. 

These two curves measure the best and the worst 
cases of the fluctuations. We found that the relative 
fluctuations for long time scales of t > 100 ms reach 
up to 3% in the worst case. They are thus one order 
of magnitude greater than the variations caused by 
intensity fluctuations. The frequency fluctuations cr^r) 
are again calculated using Eq. H48|l. This result is 
plotted together with the observed visibility in Fig. 13 
Our data points lie in between these best and worst case 
predictions. 

(3a) Heating effects. Heating processes in the trap 
can also cause significant irreversible decoherence, since 
they cause a variation of the atomic resonance frequency 
within the microwave pulse sequence. A constant heating 
rate E increases the average energy of the atoms for the 
second free precession interval [ttt, 2r7r], compared to the 
first interval [0,T,r]j by -Er^. The energy E of individual 
atoms, however, can be changed by much more than this 
average energy gain. 

To estimate the effect we have to calculate the typical 
energy change of individual atoms during the free pre- 
cession interval caused by the fluctuating forces which 
are responsible for the heating. For this purpose we 
approximate the trap as harmonic and assume the fol- 
lowing model of the heating process. Due to a random 
walk in momentum space, an initial atomic momentum 
p = \/2mE evolves into a symmetric, Gaussian momen- 
tum distribution around p with a standard deviation of 
Aprms given by the average energy gain Ett^ 



Et-, = 



2m 



(49) 



Assuming E 3> Et^^ we can linearly approximate the 
energy-momentum relationship at E. In this approxima- 
tion the momentum distribution is therefore equivalent 
to a Gaussian distribution of the energies with 



ASnns = 2\ EEt^ 



(50) 



According to Eq. (|14|l the corresponding standard devi- 
ation (Theat,_E of the dctunings AS is 



a-heat,E(T,r) = -yEEr^r, 



(51) 



depending on the initial energy E. 

We now integrate the distribution of the detunings 



Pe{AS) 



over the initial energy E weighted by the n- 
dimensional thermal energy distribution p^"-* (E) cx 
E'^^^expi-E/kBT): 



p(")(A(5) = / pe{ 
Jo 



{AS)p^''HE)dE. 



(53) 



Finally we obtain the rms detuning fluctuations (Thcat 
from the resulting distribution of AS as 



(n) 
heat 



A(S2p(")(A(5)dA<5. 



(54) 



Evaluation of (T^"^^ for the experimentally relevant time 
scale — T2/2 yields 



heat 



-ET^ksT. 



(55) 



Heating effects in our trap have been investigated in 
detail in Ref. upper limit for the heating rate 

oi E = 2 ■ 10~^ mK/s is obtained from the trap lifetime 
of 50 s (for Uq — 1.0 mK). When we linearly scale E 
to our trap depths of Uq = 1.0 mK, Uq = 0.1 mK, and 
Uq = 0.04 mK, and assume temperatures of T = 0.1 mK, 
T = 0.06 mK, and T = 0.02 mK, we obtain fluctuation 
amplitudes for the 3D-case (n = 3) of <jI^^^^ — 5.3 Hz, 

'''heat = ^-^ ^^'^ ^heat ~ respectively. We 

stress however, that these values for CThoat are upper 
limits since we did not measure heating rates E for the 
trap depths we used. The actual values of E and the 
resulting values for Chcat could be orders of magnitude 
smaller than the upper limits inferred from the life 
time because the heating rate strongly depends on the 
oscillation frequencies and the details of the laser noise 
spectrum. 

(3b) Photon recoil. Our model of the heating process 
also gives an estimate of the dephasing due to photon 
recoil. If we had one photon scattered per time interval 

giving two recoils, we would obtain a heating rate 



E : 



m 

Inserting into Eq. (|55|l {n — 3) yields 



iph 

heat 



r]k 



BksT 



Scattering of nph photons would yield 
a■ph(?^ph) = 



(56) 



(57) 



(58) 



Given a scattering rate Fg, the number of scattered 
photons obeys a Poissonian distribution. Since for our 
parameters, the probability of scattering more than one 
photon is negligible, we obtain 
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We use the temperatures of the previous paragraph and 
the photon scattering rates (see below) of Fg — 10.6 s^^, 
Fg = 1.06 s~i and F^ = 0.41 s^^. With = T^/2 
we obtain o-p\J(Fl^ — 4.5 Hz, crph(T2) = l-^ Hz and 
crph(T'2) = 1.4 Hz, respectively. 

H) Fluctuating magnetic fields. Using a fluxgate mag- 
netometer we measured a peak-to-peak value of the mag- 
netic field fluctuations of AB = 0.13 /zT, dominated by 
components at = 50 Hz. The resulting frequency shift 
on the |F = 4, mi? = 0) \ F = 3,mF = 0) transition 
is: 



BqAB, 



(60) 



where Bo = 97.9 /iT is the offset field and A(jJo^o/2vr = 
43 mHz/(/xT)^ is the quadratic Zeeman shift. For our 
case, we obtain Aw = 1.1 Hz. 

The effect of the magnetic fluctuations depends on 
the time interval between the microwave pulses. If this 
time is large compared to l/v, all fluctuations cancel 
except for those of the last oscillation period. As a 
consequence, the effect on the detuning fluctuations a 
also decreases. We calculate this effect by computing 
the Allan deviation o'a(t) of a 50 Hz sine signal. The 
detuning fluctuations then read (Thir) = ^/2Auj aA{T). 
The resulting (ibiT^), shown in Table ITTTl is too small 
to account for the decay of the spin echo amplitude. 

(5) Fluctuation of microwave power and pulse dura- 
tions. The application of two 7r/2-pulses and one tt- 
pulse results in u;ocho(2T7r) = —1. Any fluctuations of 
the amplitude (AOn/riR) or pulse duration (Ar/r) re- 
sult in variations of the amplitude of the spin echo signal, 
i.e. Wechoi'^TT^) = — cos A0 according to: 



\ 27r 



At 

T 



(61) 



With Ar/r < 10"^ (measured) and AQn/flR < IQ-^ 
(speciflcations of the synthesizer) we obtain A(j)/2TT < 
10~^, which is too small to be observed. Moreover, this 
effect neither depends on the dipole trap depth nor on 
the time delay between the microwave pulses. 

The timing of the microwave pulses would be affected 
by a clock inaccuracy of the D/A-board of the computer 
control system which triggers the microwave pulses. 
Its specified accuracy Ar/r = 10~^, results in a phase 
fluctuation S't^^At/t < 0.01 for all parameters 6' and 
T-rr used in our experiment. Thus, the fluctuations of 
microwave power, pulse duration, and timing do not 
account for the observed reduction of the spin echo 
visibility. 



Raman processes are strongly reduced due to a destruc- 
tive interference of the transition amplitudes. Thus, the 
spin relaxation rate is much larger than the spontaneous 
scattering rate. This effect was flrst observed on optically 
trapped Rubidium atoms in the group of D. Heinz en l22l 
and was also verified in experiments in our group [23 • 

The corresponding transition rate is calculated by 
means of the Kramers-Heisenberg formula 24] , which is 
a result from second order perturbation theory. We ob- 
tain for the rate of spontaneous transitions, Fg, from the 
ground state \F,m) to the ground state | F" , m" ) : 



4hd^ 



,(1/2) 



,(3/2) 



Al/2 



^3/2 



(62) 



where Aj/ = o^l — t^j' is the detuning of the dipole 
trap laser from the ^Pj/-statc and d — (4,4|/i+i |5,5), 
with the dipole operator for Am = +1 transitions. 
The transition amplitudes a''' ' are obtained by sum- 
ming over all possible intermediate states \F' ,m') of 
the relevant ^Pj' manifold For Rayleigh scatter- 

ing processes, which do not change the hyperfine state 
{F,M = F",M"), the amplitudes add up, a<'/"' = 
2a*^^^'. However, for state changing Raman processes 
{F, M 7^ F",M"), the two transition amplitudes are 
equal but have opposite sign, a'^^^' = — a'^^^'. Then 
the two terms in Eq. H62|) almost cancel in the case of 
far detuning, A1/2 ^^3/2- As a result the sponta- 
neous Raman scattering rate scales as 1/A^ whereas the 
Rayleigh scattering rate scales as 1/A^. The suppression 
factor can be expressed using the fine structure splitting 



A. 



^3/2 



Al 



Raman 



/2 
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with P 



3Ai 



/2 



(63) 



For the case of cesium, we obtain a suppression factor 
of /3 = 0.011. The Rayleigh scattering rate for an atom 
trapped in a potential of Uq = 1.0 mK is Fg — 11 s~^. 
Then, the corresponding spontaneous Raman scattering 
rate is Fpiaman = 0.12 and the population decay time 
Ti = rRaman ~ ^- Sincc in most of our experiments, 
the trap depth is significantly smaller, Ti will be even 
larger. As a consequence, we neglect the population 
decay due to spontaneous scattering. Note that the 
experiments of Refs. |2^ |23| were only sensitive to 
changes of the hyperfine P-state, since the atoms were 
in a mixture of mi?-sublevels. However, the theoretical 
treatment above predicts similarly long relaxation times 
for any particular mF-sublevel, which is consistent with 
our observations. 



(6) Spin relaxation due to light scattering. The pop- 
ulation decay time, Ti, is governed by the scattering of 
photons from the dipole trap laser, which couples the two 
hyperfine ground states via a two-photon Raman transi- 
tion. In our case, the hyperfine changing spontaneous 



C. Conclusions 

We have developed an analytical model which treats 
the various decay mechanisms of the hyperfine coherence 
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of trapped cesium atoms independently. This is justified 
by the very different time scales of the decay mechanisms 
{T2 ^ T2' <C Ti). Our model reproduces the observed 
shapes of Ramsey and spin echo signals, whose envelopes 
are the Fourier transform of the energy distribution of 
the atoms in the trap. 

The irreversible decoherence rates manifest themselves 
in the decay of the spin echo visibility and are caused 
by fluctuations of the atomic resonance frequency in 
between the microwave pulses. In the above analysis 
we have investigated various dephasing mechanisms 
and characterized them by the corresponding amplitude 
of the detuning fluctuations which are summarized 
in Table IIIII We find that a major mechanism of 
irreversible dephasing is the pointing instability of 
the dipole trap laser beams resulting in fluctuations 
of the trap depth and thus the differential light shift. 
Significant decoherence is also caused in the shallow 
dipole trap by heating due to photon scattering. Heating 
due to technical origin, such as fluctuations of the depth 
and the position of the trap, cannot be excluded as an 
additional source of decoherence. 

Compared to our experiment, significantly longer co- 



herence times (T2* = 4 s) were observed by N. Davidson 
and S. Chu in blue detuned traps in which the atoms 
are trapped at the minimum of electric fields 0. In 
Ref. 0, T2 = 15 ms obtained with sodium atoms in a 
Nd:YAG dipole trap {Uq = 0.4 mK) was reported, which 
is comparable to our observation. In other experiments, 
the inhomogeneous broadening has been reduced by the 
addition of a weak light field, spatially overlapped with 
the trapping laser field and whose frequency is tuned 
in between the two hyperfine levels Of course, 

cooling the atoms to the lowest vibrational level by using 
e.g. Raman sideband cooling techniques [26L l27| . would 
also reduce inhomogeneous broadening. The magnetic 
field fluctuations could possibly be largely suppressed by 
triggering the experiment to the 50 Hz of the power line. 
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